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The Boltzmann equation for excess Cooper pairs above T c is derived in the framework of the time- 
dependent Ginzburg-Landau (TDGL) theory using Langevin's approach of the stochastic differential 
equation. The Newton dynamic equation for the momentum-dependent drift velocity is obtained 
and the effective drag force is determined by the energy dependent life time of the metastable Cooper 
pairs. The Newton equation gives just the Drude mobility for the fixed momentum of Cooper pairs. 
It is shown that the comparison with the well-known result for Aslamazov-Larkin paraconductivity 
and BCS treatment of the excess Hall effect can give the final determination of all the coefficients of 
TDGL theory. As a result the intuitive arguments used for an interpretation of the experimental data 
for fluctuation kinetics are successively introduced. The presented simple picture of the degenerated 
Bose gas in r-approximation near the Bose-Einstein condensation temperature can be used for 
analysis of fluctuation conductivity for the cases of high frequency and external magnetic field for 
layered and bulk superconductors. The work of the Boltzmann equation is illustrated by frequency- 
dependent Aslamazov-Larkin conductivity in nanowires, in the two-dimensional case and in the case 
of strong electric field where the TDGL equation is solved directly. There are also derived explicit 
formulas for the current in the case of arbitrary time dependence of electric field up to THz range, 
the distribution of fluctuation Cooper pairs for nonparabolic dispersion, the influence of the energy 
cut-off and the self-consistent equation for the reduced temperature. The general theory is illustrated 
by formulas for fluctuation conductivity in nanowires and nanostructured superconductors. 
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I. INTRODUCTION 

For all high temperature superconductors the fluctua- 
tion phenomena can be observed and their investigation 
takes a significant part of the complete understanding of 
these materials; for a contemporary review on the fluc- 
tuation phenomena in superconductors see the review by 
Larkin and Varlamo\i The Ginzburg-Landau (GL) ap- 
proach of the order parameter is an adequate tool to in- 
vestigate the low-frequency behavior of fluctuations near 
to Tr\ for a review of the Gaussian GL fluctuations see 
Ref. |3- A lot of important papers on the fluctuation phe- 
nomena in superconductors and related topics have not 
been cited in these reviews, see for example Ref. 0. We 
wish to point out that the GL approach is the standard 
tool for the investigation of magnetic field penetration 
in superconductors^ and even non-Gaussian approach to 
critical fluctuations £ 

Amidst all kinetic phenomena the fluctuation conduc- 
tivity created by the metastable in the normal state 
Cooper pairs is probably best investigated. The Boltz- 
mann equation is a standard tool for investigation of ki- 
netic phenomena and the purpose of the present paper is 
to derive the Boltzmann equation for fluctuation Cooper 
pairs and to illustrate its work on the example of the fluc- 
tuation conductivity; a shortened version of the present 
research was presented in preliminary communications^ 
We rederive the frequency dependence of the Aslamazov- 
Larkin conductivity, fluctuation Hall effect at weak mag- 
netic fields, and magnetoconductivity. We analyze the 



experimental data for indium oxide films and find signif- 
icant deviation from the BCS weak coupling prediction. 
We are coming to the conclusion that a systematic in- 
vestigation of lifetime of fluctuation Cooper pairs will 
give important information for our understanding of the 
physics of superconductivity 

II. FROM TDGL EQUATION VIA BOLTZMANN 
EQUATION TO NEWTON EQUATION 

Our starting point is the time-dependent Ginzburg 
Landau (TDGL) equation for the superconducting order 
parameter derived in the classical paper by Abrahams 
and TsunetO)i see also Ref. and references cited in the 
review by Larkin and Varlamo\ii 

2m J ^ + a^> + b\^\H^-h 1 (D t ^-0, (1) 

where m* and |e*| = 2|e| are the mass and charge of the 
Cooper pairs, parameter 7 describes the dissipation, and 
£(r, t) is the external noise in TDGL equation. Here 

- ihD r = -ih\7 - e*A/c, 
iHD t = ihd t — e* 'ip, 

are the operators of kinetic momentum and energy, A is 
a vector-potential, and tp is the potential. 

Close to the critical temperature a(T) « (T—T c )ao/T c , 
and b const, where ao = ft 2 /2m*£ 2 (0), and £(0) is 
coherence length. 
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The correlations of the white noise (Q = 0, 

(C*(ri,ti)C(ra,*a))=r5(ti-t a )5(ri-ra), (2) 

are parameterized by fluctuation parameter T. The BCS 
theory gives 



7i 



7T Oo 

8 T ' 



and that is why we parameterize 7 = 7 BCS r re ;, by the 
dimensionless parameter r re i ~ 1, which describes the 
relative life-time of fluctuation Cooper pairs. 

The most simple is the case of free particle, which 
means A = 0, tp = 0, 6|'I , | 2 ss 0. Introducing the Fourier 
transformation 



ip-r/H 



ip-r/H 



c('»*) = E— at^*)' 



where 



and 



v/V 
(2tt^) £ 



(3) 
(4) 



(^(*i)^(*2))=r« M *(t 1 -t a ), 
we obtain TDGL equation in momentum representation 
(e p + a)V> P = -hj(d t ip p - ( p ). (5) 
The solution of this reads 

VVW = e-* /2Tp (^j e*'/^^ Cp( i')^ + ^(q)) , ( 6 ) 



where 



^0 



2(e p + a(T)) ' 



2m* 



(7) 



are momentum-dependent lifetime and kinetic energy of 
fluctuation Cooper pairs. The number of particles for 
every momentum can be found by noise averaging 

n p = (r p (W P (t)) = np(0)e-*/ T » + (1 - e^Jfip, (8) 

where n p (0) = \ip p (0)\ 2 is the initial number. The time 
differentiation of this solution gives the well-known Boltz- 
mann equation 



^-n p (t) = - — (n p (t) - n p ), 
at t p 



(9) 



which can be considered in this physical situation as a 
consequence of the TDGL equation. The quantity 

h p — n p (t = 00) = IYp 



gives the equilibrium number of particles. The fluctua- 
tion parameter T is related to dissipation parameter 7 by 
the fluctuation-dissipation theorem, which here takes the 
form 



r _2T _ ?i p _ T 



o-qTo 



(10) 



where 



T 



P e p + a{T) 

is the Rayleigh-Jeans distribution. 

Let us now analyze the influence of a weak electric field 
in the BoltzmannS equation 



dfUp + e*E • dpn p 



1 



(rip - n p ). 



(11) 



Quantum mechanics was born during Halle confer- 
ence in 1891, when exposed to the ignorant criticism of 
both statistical methods and atomic physics, Boltzmann 
suddenly made a remark: "I sec no reason why ene rgy 
shouldn't also be regarded as divided atomically" Ref . Ilfl 
Later on applying the Boltzmann method to the problem 
of black body radiation Planck found that the constant 
appearing in the photon spectrum is just the volume of 
the Boltzmann cells in the phase space. Due to this rea- 
son, Planck called the quantity 2nh after Boltzmann - 
Boltzmann constant. 

For the solution we search in the form 

n J ,(i)=»(p,*)wn(p-m*V(p,t)) ) (12) 
and we obtain the Newton equation 



m*d t V p (t) = e*E V p (t) 



(13) 



for the field of drift velocity in momentum space. The 
general formula for the current gives 



V 



where is the D-dimcnsional volume density of the 
fluctuation Cooper pairs. Substitution here of the shifted 
equilibrium distribution gives the well-known formula for 
the conductivity tensor— 



an = e 



d D p v p ® v p / d n p 
(2irh) D 1/t p - iu \~d£p~ 



(15) 



where v p = d p e p = p/m* is the Cooper pairs' velocity. 
This is only a small fraction of the total conductivity 

T T — T ' 

^(T) = ° N (T)+M e ), 6 = In — » - c , cr fl < < T N . 

(16) 
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For thin superconducting films D = 2 substituting 



integral 



dp x dpy d(irp 2 ) m* dn T 



(17) 



{2ttK) 2 {2ttK) 2 2ttH 2 p ' de (e + a) 2 ' 
we obtain the classical result by Aslamazov and Larking 



16h Trel T-T r 



irh 2 



where 



r(e) =r(p = 0,e) = 



Trfi- r rc i 



Tie) 



TO 



(18) 



(19) 



is the lifetime for Cooper pairs with zero momentum. 

For the two-dimensional (2D) case conductivity is just 
the inverse resistance a^ 2D \T) = i?^ 1 (T). For con- 
ventional disordered superconductors normal conductiv- 
ity can be approximated by residual conductivity far 
above T c , for example, T = 3T C . In this approximation 
Aslamazov-Larkin conductivity can be rewritten in a con- 
venient for experimental data processing form 



R n (T) R a (3T c 



I6h (T 



e 2 T re l \T, 



(20) 



Performing the linear regression fit of the data pre- 
sented in Ref. ^| we have obtained that for indium oxide 
films T rc i = 1.15. This significant 15 % deviation from 
the weak coupling BCS value is created by strong cou- 
pling effects. Wc conclude that analogous systematic in- 
vestigations for thin films would be very helpful for our 
understanding of the dynamics of the order parameter 
in superconductors. Decreasing the lifetime and r rc i by 
depairing impurities or disorder for anisotropic gap su- 
perconductors definitely deserves a great attention. 



III. FLUCTUATION CONDUCTIVITY IN 
DIFFERENT PHYSICAL CONDITION 

A. High frequency conductivity 

For diagonal components of conductivity taking into 
account that Trl = D, from the general formula Eq. (|15fl . 
we obtain 



-,*2 



D 



d D p 



{2ttK) d 1/t p - iuj 



de T , 



(21) 



It is convenient to introduce a dimensionless frequency 
z = ujT(e). In order to derive the dimensionless com- 
plex conductivity ^(lo) we need to solve the elementary 



x 2 (x + y) 



dx 



1 + - ln(l + i/)-l 

y 



(22) 
(23) 



1 H ln(l - iz) - 1 



where x = p 2 /2m*a(e) + 1 is the kinetic energy of Cooper 
pairs taken into account from the "chemical potential" in 
o(e) units, and y — —iz = —iujT(e) is the dimensionless 
Matsubara frequency C M = — The integral Eq. 12 I'll is 
solved considering the Matsubara frequency y to be a real 
variable. Then we can make the analytical continuation 
to real frequencies substituting y = —iz in the result 
Eq. (|23|l . This method is very popular in the quantum 
field theory, but works effectively for classical problems 
as well. In such a way we obtain 



^( Z ) = 72 



-V 



z 
2 

TT 

2 r 

2 



zarctan(z) — — ln(l 



3^2 (x) 



dx, 



z 



arctan(z) — z + — ln(l + z 2 ) 



2 ^ 
— V 

7T 



dx. 



x* — z^ 

Then the frequency-dependent conductivity reads 
o-2d( £ ^) = o- AL (e)^(^T(e)). 



(24) 



(25) 



(26) 



The integral Eq. JSTJ can be solved for arbitrary dimen- 
sion 



cr D (e,cj) = (T D (e)<r D (z), 
cf. the paper by Dorsej^ 



(27) 



i \ , r(2-^/2) e 2 2 _ c Tr(e) 
^ D (e)-4 (4?r)D/2 -K(e)] — , (28) 

where £(e) = £(0)/y/e is the temperature-dependent co- 
herence length. The conductivity in this case has the 
form 



^i,d( z ) 



?2,d(«) 



[1 



D(D-2)z 2 
-(l + z 2 ) D / 4 



D{D-2)z 2 



D 

I — arctan z 
D 



+(1 + 



2ND/4 



D > 
I — arctan z 



(29) 
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B. Hall effect 



D. Strong electric fields 



The fluctuation Hall conductivity also can be derived 
in the framework of the Boltzmann kinetic equation. We 
have to take into account a small imaginary part a of 7 
parameter in the TDGL equation, i.e., 7^7 + ia, and 
a -C 7. The solution^ of the kinetic equation gives 
y 

Vxyit) = -w c r(e)(JAL(e) oc r 2 (e), (30) 

where uj c = e*B/m*c is the "cyclotron" frequency and 

1 



Z = -Im- 



7' 



a 

— « L 

T 



(31) 



This result agrees with microscopic calculations^ Due to 
the small value of the parameter a, fluctuation Hall effect 
is difficult to observe. With fitting of a and m* from the 
experimental data finishes the complete determination of 
parameters of TDGL theory 



C. Magnetoconductivity 

It is interesting to mention that the classical formula 
for the conductivity Eq. I|15|) correctly works even for 
strong magnetic fields. We only have to substitute the 
momentum integration with summation on discrete Lan- 
dau levels, taking into account the density of Landau 
magnetic subbands 



where 



= hi 



h 



huo c 

2a 



- ) =ao(2n+l)h, 



B c2 (0) 

is the dimcnsionlcss magnetic field and 

d 



B c2 (0) 



Trj^B c2 {T)\ Tc 



(32) 



(33) 



(34) 



is the linear extrapolation. In the numerator of Eq. I|15|) 
we have to substitute the classical velocity with the os- 
cillator matrix elements of the momentum. Analogously 
for the energy-dependent lifetime we have to average on 
neighboring levels. Due to the triviality of the oscilla- 
tor problem these substitutions can be performed in only 
one way, and Aslamazov-Larkin conductivity Eq. i|18|) 
is substituted by the magnetoconductivity of Abrahams, 
Prange and Stephen^ (APS) 



e 2 Tr n 1 



> CT, 



, s (e,h) 



e 2 Tr 



i.e., 1/e has to be substituted by APS function 



f(e,h), (35) 



f(e,h) = 



h 2 



cf 



c 

2h 



-ef 1 



e 

2h 



(36) 



This two-dimensional result can be easily generalized for 
layered and bulk superconductors using the layering op- 
erator introduced in Ref. 



Using the optical gauge 

ip = 0, A = -fE, 
the TDGL equation Eq. Q reads 

d u %{u) = -i [(q + fuf + e] il> g (u) + C,(u), (37) 

where we are introducing dimcnsionlcss variables for the 
u = i/Vo time, q = p£(0)/h momentum, / = e* Etq^(0) / h 
electric field, and ( q (u) = r ( p (t) noise. We have a lin- 
ear ordinary differential equation which can be solved for 
arbitrary f(u), i.e., for arbitrary time dependence of the 
electric field. For constant electric field the TDGL equa- 
tion has the solution 



ip q (u) 



exp 



[(q + fu 2 ) 



du-. 



x exp 



X ( q (ui)dlli + Vg(0)} 



[(q + fusf + e] du 3 



(38) 



In order to obtain the static [t 3> t(c)] momentum dis 
tribution we have to perform the noise averaging 

n k = lim (\ip q+fu (u)\ 2 ) 



exp 



T 

a-o Jo 

where u\ = u — v and 



1 



2„,3 



(k 2 + ( )v + fkv 2 - -fv 



(39) 



dv, 



k = q + fu = (p- e*A) 



e(o) 



is the dimensionless kinetic momentum. This distribu- 
tion can be directly derived from Boltzmann equation 
Eq. lfTT|) for fluctuation Cooper pairs w& In Ref. Ha it was 
demonstrated that substitution of the momentum distri- 
bution Eq. (|39[l in the formula for the current density 
Eq. 114f> gives the result which agrees with the formula 
by Dorseyi^; cf. also the paper by Gor'ko\ii£ 



3 {Ex) 



16fi [ZvrVaf (0)] 



D-2 



exp 1 



gu 3 ) 



,(D-2)/2 



du, 



(40) 



where 



9 = Tr f 



e*E x £(0)T _ it eE x £(0) 



T re l. 



Differentiating the upper expression we obtain differen- 
tial conductivity 



Cdiff 



dj(E x ) 
dE x 



e 2 r re i 



1 



\Qh{2^{Q)) D - 2 
- 2gu 3 



,(D- 2 )/ 2 



exp(— eu — gu ) du. (41) 
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Applying a voltage U(t) = U DC + U AC cos ut to the 
nanowire, the differential conductivity can be easily de- 
termined measuring the AC component for the current if 
U AC <C U DC . Cooling the sample the differential conduc- 
tivity will decrease, then at some temperature it will be 
annulated and what will happen at further cooling is an 
interesting experimental question. 



IV. CURRENT FUNCTIONAL: 
SELF-CONSISTENT APPROXIMATION AND 
ENERGY CUT-OFF 

The self-consistent approximation for the reduced crit- 
ical tempcrature^iia in the one-dimensional (ID) case 
reads 

£ro„ = In — + — n 1D = In ^ + e 1G M (e rcn , /), (42) 

where n 1D is the bulk density of the fluctuation Cooper 
pairs when we have ID fluctuations in a wire with cross 
section S <C £ 2 (e) and 



M oA 2 (0K(0)e 2 r c 



80FAC£(O)S' 



(43) 



where A(e) = A(0)/y— e is the temperature-dependent 
penetration depth and AC is the jump of the specific 
heat at T c per unit volume. For numerical calculations 
the function 



A/i(e,/) 



exp(— ev — gv ) 



dv 



has to be programmed as 

/>oo 

M(e, /) = 2 / cxp(-gz 6 - ez 2 )dz 
Jo 



(44) 



(45) 



Analogously for the thin superconducting film with thick- 
ness df <C £(e) the equation for reduced temperature at 
zero electric field takes the form£ 

£ren = m ^ + — n 2D = In + 6 2G AA 2 (e rc „), (46) 
J o a o 1 o 

where n 2D is the volume density of the fluctuation Cooper 
pairs having 2D fluctuations, 

^-jsz^-HK^)' (47) 

is the 2D Ginzburg number and 

'c + e 



Af 2 (e) = In 



(48) 



As simplest possible application of these results we 
have to mention nanostructured superconductors, e.g., 
nanowires^, similar to those used for long time for in- 
vestigation of phase slip centers^ of the superconducting 



phase. We are pointing out that paraconductivity is a 
property of the normal phase. 

For general (nonparabolic) dispersion we can derive 
from the TDGL equation the formula for the distribution 
of fluctuation Cooper pairs 



n k (u) = — [ expi- / 
a-o Jo I Jui 

+ n k (0)exp 



+ e 



du-2 > du\ 



a 

e{k{u 2 )) 
a 



where the dimensionless kinetic momentum and the vec- 
tor potential are 



du 2 \,(49) 



k(u) — q + A(u), A(u) 



A(t). 



In the case of parabolic dispersion and arbitrary time 
dependence of the electric field we can write Eq. I|14|l for 
the current functional in the form 



— f F A [m]dui+n fe (0)F A [0] 
ao Jo 



m*^(0) 

where for brevity we introduce the functionals 
A(u) B A [m 



F A [ui] 



y/u — U\ (u — Ul) 3 / 2 

{Ba[uA) 2 



x exp 



u — U\ 



- G{ui) 



(50) 

(51) 
(52) 



and 



B A [ui] 



A(u 2 )du 2 , 



G A [ 



e(u-iti)+ / (A(u 2 )) du 2 



Local GL theory with an energy cut-off for the kinetic 
energy e p < a c [or for the kinetic energy taken into ac- 
count from the chemical potential e p +a(T) < aoc] is very 
often used for describing the fluctuation phenomena™ 



a k 2 , \k\ < A 
a c, \k I > A 



(53) 



where the dimensionless constant A = y/c ~ 1. 
For YBa2Cu307_5 the recent investigations^ , of high- 
frequency fluctuation conductivity determined A sw 0.5. 
See also the recent work by Puica and Lang, Ref. 0, 
where the effects of energy cut-off and self-consistent in- 
teraction are treated in great details. Then the functional 
which participates in the formula for the current Eq. I|50|) 
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takes the form 



F A [ui] 



smh.(2AB A [u 1 ]) 

y/u — Ui 



x exp 



— (u — Ml) C + 



[U - Ml) 



(#aM) 2 

(it — Ml) 2 
U — Ml 



(54) 



r A-A(u) 

exp 

-A-A(u) 



(m - Mi) g + 



M — Ml 



These rather complicated formulas are necessary for in- 
vestigation of paraconductivity in the THz range. In the 
next section we will give an illustration for the important 
one dimensional case. 



V. FLUCTUATION CONDUCTIVITY IN 
NANOWIRES 




FIG. 1: Age distribution for fluctuation Cooper pairs above T c 
according to electric field f e . At absciss is the momentum k e and 
at ordinate is the function J-"_|_(fc e ,/ e ) from Eq. 1581 . The curves 
illustrate different values of / e . The number of particles with zero 
kinetic momentum in weak electric fields is much more than in 
strong fields. 



The recent development of the technology of the per- 
formance of nanowires made it possible and even indis- 
pensable for the investigation of fluctuation conductivity. 
In this section we will analyze in detail the general results 
in ID case. 

The integrants in the momentum distribution Eq. 139|) 
is actually age distribution for fluctuation Cooper pairs 



f{v; k,e,f) = exp 



(k 2 + e) v + fkv 2 - -/V 
3 



(55) 

the variable v is the age in units r(e). Time integration 
returns us to the momentum distribution, which using 
the dimensional variables 



reads 



where 



11/2 > 



n(k; ej) 



f 



13/2 



F±(ke,fe), 



(56) 



(57) 



F±(k e , / e ) = / exp 
Jo 



-(ki±l)x + f e k e x*-^-f*x 3 



and x = \e\u. 

For the case e = 0, T = T c or when / — > oo, i.e. 



dx 
(58) 



U = e*E x i{e)r{e)/h » 1, f(e) = e(0)/|e| 1/2 (59) 
we obtain 




FIG. 2: Age distribution for fluctuation Cooper pairs under T c 
depending on electric field / £ . At absciss is the momentum fc e and 
at ordinate is the function J-"_(fc e ,/ E ) from Eq. 1581 . The curves 
illustrate different values of / e . Decreasing the electric field the 
number of particles increases. The maximum in strong electric 
fields is right-shifted for particles with positive momentum. 



where in 



exp 



-kjy + k f y 2 - -y 3 



dy (61) 



we use the transformation y = f 2 ^ 3 v. 

For ID case using Eq. (|4(J|I . we express the current 



i(kj) 



riT 



(60) 



~~8h~ 



(62) 
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FIG. 3: Momentum distribution for fluctuation Cooper pairs in a 
strong electric field at T c . At absciss is the scaling momentum kf 
and at ordinate is the universal function Ta(kf) from Eq. 1611 . 
which describes particle distribution. 



FIG. 5: Distribution of dimensionlcss function <;_|_(/ e ) according 
to Eq. 1651 above critical temperature T c . At absciss is / e . The 
maximum of fluctuational current is in zero electric field near T c . 
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FIG. 4: At ordinate is function J(e) participating in Eq. 1631 . 
which describes the strength of the fluctuational current, and at 
absciss is the dimensionlcss temperature e. The curves show the 
dimensionlcss electric field /. Coming nearer to critical temperature 
T c the fluctuational current grows up rapidly in very weak electric 
fields. 



where 



J cxp(-ew - gv 3 )y/vdv, 9=^- (63) 



The fluctuational current for Cooper pairs above and 
under the critical temperature is 



3 2 T rc i£(0)£ a 
16ft|e| 3 / 2 



?±(/e), 



(64) 



FIG. 6: Asimptotic J* <;-|_(/ e ) from Eq. 1661 is going to constant 
when / £ — > oo. Above T c in very strong electric field fluctuational 
current goes to constant. 



where 



?±(/e) 



2 f°° 

-= / exp(=F^ £ - g e v e )Jv^dv e 
/7r Jo 



(65) 



is the dimensionlcss function, which depends on the 
strength of the electric field. For convenience we use 
g t = / e 2 /f2 and v e = v\e\. We wish to point out the nor- 
malization and the asymptotics 



?+(0) = l, ?±UWoc) 



V3/e 



(66) 



Using the strong field asymptotic f e ^> I, we obtain 
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FIG. 7: Near f e -* 0, asimptotic ^|^ f+ (/ £ ) from Eq. JHEJ, which 
describes fluctuational current above T c , is not sharply-outlined, 
but has a smooth decrease. The fluctuational current is small at 
high temperatures. 
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FIG. 8: Function A/i(e,/) from Eq. 1711 is at ordinate and e - at 
absciss. The curves illustrate the electric field / near the phase 
transition. Very close to the critical temperature T c the ID di- 
mensional density of Cooper pairs has great fluctuations in strong 
electric field and small ones in weak electric field. 



the fluctuational current at T c 



j(/e -> OO) 



_2_eJ; 
V3 ti ' 



(67) 



We wish to point out that j(f e — > oo)/T c is universal 
and contains only fundamental physical constants. All 
material constants like £(0), A(0) or the cross section of 
the nanowire S are cancelled. For experimental data pro- 
cessing it is necessary to perform linear regression of the 
IV curve 



Jtot 



U 



Rn(T c ) 



+ sign((7) x 24.22 nA T C [K], (68) 



where the second term is the universal fluctuational cur- 
rent at strong electric fields around T c . In order to avoid 
the thermoelectric effect coming from different materi- 
als forming the contacts to the nanowire one can an- 
alyze the current harmonics predicted by Eq. 168fl for 
U(t) = U smiot at T = T c 




FIG. 9: At ordinate is A/i(e,/) from Eq. 1711 and at absciss is e. 
The ID dimensional density decreases with increase the tempera- 
ture. The curves illustrate the electric field /. 



Jtot(i) 



Ua 



Rn{T c 



sin tot H — j c 

7T 



2 = 1 



sin[(2Z + l)tot] 
27+1 



(69) 



This current response is a good approximation even 
slightly above T c for voltage amplitudes 



8T e L 
it e £(0)' 



|3/2 



(70) 



where L is the length of the nanowire. In such a way the 
investigation of fluctuation current in nanowires can be 
used as a high accuracy test for applicability of the TDGL 
equation for nanostructured superconductors. For fur- 
ther references related to harmonic generations close to 
T c see Ref . EIH 



Using Eq. 114fl and Eq. (|39|l we derive ID dimensional 
density 



niD 



riT 



2^(0)5 



M(e,/), 



(71) 



where A/"i(e, /) is previously defined in Eq. 1)44(1 . when we 
considered without a derivation the self-consistent equa- 
tion for the reduced temperature. The analysis of tem- 
perature and electric field dependence is reduced to two 
functions of one variable / e above and below the T c 



TlT 



2£(0)S 



N±(f e ), 



(72) 
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FIG. 10: Illustration of M+(f e ) from Eq. J73J above T c . At ab- 
sciss is temperature-electric dependent variable f e . The density of 
Cooper pairs in nanowire is maximum in zero electric field. 
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VI. DISCUSSION AND CONCLUSION 

Solving in parallel the TDGL equation and the Boltz- 
mann equation we obtained coinciding results: not only 
for the linear case of Aslamazov-Larkin conductivity, 
but for the cases of strong electric fields, arbitrary time 
dependence of the electric field, nonparabolic momen- 
tum dependence of energy of Cooper pairs, energy cut- 
off, self-consistent equation for the renormalizcd reduced 
temperature, frequency dependence of the fluctuation 
conductivity etc. The number of fluctuation Cooper 
pairs which participates in the Boltzmann equation and 
the formulas for the current is actually the diagonal 
element of the order parameter correlators^ rifc(t) = 
C [pkin = P — e*A(t);t, t] . One can also easily check that 
the entropy of fluctuation Cooper pairs r\ is increasing 
with the time drj/dt > 0; the capital rj in the 77-theorem 
by Boltzmann is often spelled as Latin H. Our self- 
consistent formula for the fluctuation conductivity of a 
superconducting nanowire can be directly used for the 
experimental data processing. In such a way we con- 
clude that Boltzmann equation for fluctuation Cooper 
pairs reproduces the known results of the fluctuation the- 
ory in the normal phase and it is a adequate tool to pre- 
dict new phenomena related to mctastable Cooper pairs, 
like negative differential conductivity in the fluctuation 
regime predicted in Ref . ^| and strong electric field effect 
in nanostructured superconductors where the heating ef- 
fects are reduced. Our universal result Eq. (jfiSJI for a 
fluctuational current in a nanowire under strong electric 
field shows that Boltzmann equation for the fluctuation 
Cooper pairs will become an indispensable tool for un- 
derstanding the electronic processes in nanostructured 
superconductors. 



FIG. 11: Asimptotic A/+(/ e ) from Eq. 1731 . when f € — ► crosses 
smoothly ordinate at 1. At absciss is temperature-electric depen- 
dent variable / e . 



where 

A4(/e) 



N+(f t = 0) = 1. 

(73) 



Analogously at T c , where e = and / £ > 1 we obtain 
r(l/6)n T 1 



nw(f) 



(74) 



22/335/6^(0)5 /V3' 

Choosing a typical set of parameters for the Sn nanowire: 
£(0) = 1000 A, A(0) = 1000 A, D = 500 A, S = ttL> 2 /4 
according Eq. l|43|l we obtain for ID Ginzburg number 
€\g = 4.3 x 10~ 5 . This parameter is essential for the nu- 
merical solution of the equation for renormalizcd reduced 
temperature Eq. (|4*2|) . 
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